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A YIELD MANAGEMENT METHOD AND SYSTEM 

The present invention relates to a yield management method and system. 

Yield Management Systems (YMSs) are commonly used as computer-based 
supports to facilitate decisions in several applications, such as in the 
transportation business and, in particular, in. the airline industry. The 
goal of a yield management system is that of allocating an offered 
capacity to different types of requests (categories) , in order to optimise 
a yield parameter. For example, the yield management system is used to 
maximise the revenue that can be obtained from selling a given 
transportable weight and volume of a cargo flight; the objective is to fly 
an aircraft as profitable as possible without allowing low-fare freights 
to displace high- fare freights. 

All yield management systems use a forecast method to estimate the 
requests of the different categories; a decisional system is th^n applied 
to the estimated requests for determining a suggestion that should be 
given to maximise the yield parameter. Typically, stochastic or 
probabilistic methods are employed, in order to take into account an 
uncertainty of the estimated requests (and not only their mean value) . 

A solution known in the art consists of estimating a simple 
probabilistic distribution of the variables involved in the yield , 
management system. A stochastic linear programming model is then applied 
using the simple probabilistic distribution to determine an allocation 
model that maximizes a target function defining the yield parameter. An 
example of the aforementioned method is described in Williamson, E. L. 
"Airline Network Seat Inventory Control: Methodologies and' Revenue 
Impact", Massachusetts Institute of Technology Phd Thesis - Cambridge 
(Mass) June 1992. 

A drawback of the solutions known in the art is that they become 
computationally intractable when a complexity of the problem defined by 
the yield management system increases. This is particular critical if the 
offered capacity is defined by more than one variable, such as the weight 
and the volume of a cargo flight. In these cases many simplifications, 
shortcuts and assumptions are introduced, in order to solve the problem in 
a reasonable amount of time. For example, the number of categories is 
limited or a manageable form of the probabilistic distributions is assumed 
(for example a Gauss, Poisson or uniform distribution). Moreover, the 
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known solutions discard some capacity variables (for example using only 
the weight and just checking a posteriori if volume constraints are 
satisfied) or discard the stochastic nature of the yield parameter (for 
example considering only its mean value) . 

5 

However, all the simplifications described above strongly reduce the 
accuracy of the yield management system. As a consequence, the solutions 
known in the art may lead to a relevant- economic, lost. 

10 It is an object of the present invention to overcome the above- 

mentioned drawbacks. In order to achieve this object, a method as set out 
in the first claim is proposed. 

Briefly, the present invention provides a yield management method for 

15 optimising a yield parameter resulting from assigning a capacity offered 
by a future instance of a service to each one of a plurality of different 
categories- of requests competing for the capacity, the capacity l^eing 
defined by at least one capacity variable, the method including the steps 
of: storing a set of historical profiles for each one of a plurality of 

20 previous instances of the service, the set including a historical profile 
of an historical value of each capacity variable reserved by each 
category, assigning a probability to each previous instance of the 
service, estimating a potential profile of a potential value of the 
capacity variable from each historical profile according to a 

25 corresponding current value of the capacity variable reserved by the 
category for the future instance of the service and according to a 
corresponding unconstrained demand of the capacity variable for the 
category in the previous instance of the service, defining a historical 
scenario for each previous instance of the service, the historical 

30 scenario including a final potential capacity variable from each 

corresponding potential profile, and determining an authorisation to 
allocate €ne offered capacity for each capacity variable of each category 
in the future instance of the service by applying a stochastic model to 
the historical scenarios according to the corresponding probabilities. 

35 > 

Moreover, the present invention' also provides a computer program 
application for performing the method, a program product storing the 
program application, and a corresponding yield management system. 



40 



Further features and the advantages of the solution according to the 
present invention will be made clear by the following description of a 
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preferred embodiment thereof, given purely by way of a non-restrictive 
indication, with reference to the attached figures, in which: 

Fig.l is basic block diagram of a data processing system in which the 
yield management method of the invention can be used; 

Fig. 2 shows a partial content of a working memory included in the 
system; . ** . 

Figg.3a-3d are a flow chart describing the logic of the yield 
management method implemented in the system; 

Figg.4a-4b depict different diagrams of some variables involved in 
the method. 

With reference in particular to Fig.l, there is shown a data 
processing system 100, which is used as a computerised support for 
controlling cargo flights of an airline. The system 100 is formed by a 
reservation subsystem 105 and an authorisation subsystem .110; the 
subsystems 105 and 110 are connected to each other through a communication 
channel 115. 

r 

The reservation subsystem 105 consists of a mainframe having a 
Central Processing Unit (CPU) 120. A working memory (RAM) 125 and a bulk 
memory 130 (generally consisting of one or more hard-disks) are associated 
with the CPU 120. Moreover, a telecommunication control unit (TCU) 135 is 
used to connect multiple remote terminals 140 and the communication 
channel 115 to the CPU 120. The reservation subsystem 105 hosts a 
transactional processing environment, which allows the remote terminals 
140 to respond immediately to each request entered by an operator. 

The authorisation subsystem 110 consists of a workstation having a 
RISC architecture. The authorisation subsystem 110 includes a Central 
Processing Unit (CPU) 150 and a working memory (RAM) 155 that is used 
directly by the CPU 150. The CPU 150 is associated with a bulk memory 
consisting of a hard-disk 160 and of a driver unit 165 for reading CD-ROMs 
170. Moreover, the CPU 150 is connected to an input unit (IN) 175, which 
consists for example of a keyboard and a mouse, and to an output unit 
(OUT) 180, which consists for example of a monitor and a printer. A 
network interface card (NIC) 185 is used to connect the communication 
channel 115 to the CPU 150. 



4 



The above-described architecture implements a yield management 
system. Yield management systems are a well-known application field of 
several Operation Research (OR) techniques. Basically, the goal of a yield 
management system is that of optimising one or more parameters that can be 
5 obtained from a given capacity of a generic service. A demand consisting 
of multiple types of requests (categories or classes) competes for the 
offered capacity; each category is defined by a subset of the requests 
exhibiting similar characteristics (such as a physical property, a fare, a*^ 
reservation/cancellation behaviour) . 

10 

The yield management system selects the requests by accepting only 
the more convenient ones; this is done by assigning a booking limit 
(authorisation) to each category. The offered capacity may. be allocated to 
the different categories using several techniques. In an out of nesting 

15 policy, the offered capacity is simply split among the categories 

according to their authorisations. This is a very conservative technique, 
which ensures a pre-set capacity to each category; however, it m#y lead to 
a yield loss when the requests for a category are lower than the 
respective authorisation. Conversely, in a nesting policy the same 

20 capacity may be reserved by more than one category; the categories are 

ranked according to their value, in order to assign a greater capacity to 
the most valuable ones. Particularly, in a total nesting policy each 
category may' reserve all the capacity assigned to the category itself and 
to the categories with a lower value. This is a very aggressive technique, 

25 which favours the most valuable categories to the detriment of the other 
ones; however, it provides good results only when the requests for low- 
fare freights are received before the ones for high- fare freights. 
Otherwise, in a partial nesting policy each category may reserve the 
capacity assigned to the categories with a lower value only when the 

30 capacity assigned to the category itself has been completely allocated. 

In ££e example at issue, the yield management system is used to 
maximise a revenue A that can be obtained from a given capacity that is 
offered by each cargo flight; the offered capacity is defined by two 

35 different variables consisting of a total weight P tot and of a total volume 
V tot of the freights that may be transported by the flight. The freights are 
classified in five categories Ci (with i=l~5) , for example defined as 
normal light, special light, normal heavy, special heavy, and extra; this 
number is a good compromise between the opposed requirements of having 

40 uniform characteristics within each category (leading to a high number of 
categories) and a constant behaviour of each category (leading to a low 
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number of categories). 



The yield management system employs a partial nesting policy, with 



two sets of independent authorisations. A set of nesting authorisations Ap A 
is used for the weight, with the categories that are ranked in a 
respective descending nesting order ORDp (from 1 to 5) ; a set of nesting 
authorisations Av A is used for the volume, with the categories that are 
ranked in another independent nesting order ORDv. The partial nesting 
policy is preferred in this application since the period on which each 
flight is open for reservation is short, so that the order of the requests 
is very random. Moreover, it is generally difficult to rank the freights 
according to their value in an accurate manner, because both the weight 
and the volume are to be taken into account. In addition, the weight and 
the volume are not always constraints limiting the offered capacity at the 
same time. 

As described in the following, the weight nesting authorisations Ap ± , 
the volume nesting, authorisations Av if the weight nesting order ORDp and 
the volume nesting order ORDv for a future instance of each flight are 
determined by the authorisation subsystem 110. The authorisation subsystem 
110 further supplies a minimum selling price, which consists of a couple 
of values indicating a weight unitary revenue MXp and a volume unitary 
revenue MAv for each segment of a network supported by the airline, under 
which a request should be refused because not profitable. The nesting 
authorisations Api,AVi, the nesting orders ORDp, ORDv and the minimum 
selling price MXp,MXv are transmitted to the reservation subsystem 105 and 
stored on its hard disk 130. The hard disk 130 further stores booking 
information BK about every request accepted by the reservation subsystem 
105 for the future instance of each flight, such as a date of the accepted 
request, a weight and a volume of the freights, a name of the client, an 
applied fajre (if already available) . The sum of the weight and volume of 
the accepted requests gives a weight Rpi and a volume Rv ± , respectively, 
already reserved by each i-th category. The reservation subsystem 105 
calculates a residual weight Dpi and a residual volume Dv ± that are still 
available for the i-th category according to the following formulas: 
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The operator at a remote terminal 140 enters a reservation request 
for a particular flight. The reservation request is accepted if both the 
weight and the volume of the freights do not exceed the residual weight Dp ± 
and the residual volume Dvi, respectively, of the corresponding category; 
5 otherwise, the reservation request is denied. At the same time, the 

reservation subsystem 105 adds together the minimum selling price of all 
the segments from an origin to a destination of the request, and 
multiplies the result by the : weight and- the volume of- the freights; when - 
the revenue of the request (if available) is lower than the value so 

10 obtained, the request is always refused. When the reservation request is 
accepted, it is recorded on the hard disk 130 and the respective residual 
capacities Dpi,DVi are accordingly updated. The booking information BK, the 
reserved capacities Rpi,RVi and the residual capacities Dpi,DVi for the 
future instance of each flight are sent to the authorisation subsystem 110 

15 either periodically (for example at the end of each day) or on request. 

Similar considerations apply if the data processing system^has a 
different, architecture, for example consisting of a single computer or 
with the reservation subsystem and the authorisation subsystem working in 

20 the INTERNET, if the reservation subsystem and the authorisation subsystem 
have a different structure, for example if they consist of Personal 
Computers (PCs) or they are implemented with a local area network, if the 
freights are classified in a different number of categories, if other 
categories are envisaged, if a different nesting policy is employed (for 

25 example a mixed partial one), and so on. 

Considering now Fig. 2, there is shown a partial content of the 
working memory 155 of the authorisation subsystem; the information 
(programs and data) is typically stored on the hard-disk ahd loaded (at 
30 least partially) into the working memory 155 when the programs are 

running, together with an operating system and other application programs 
(not showrTin the figure). The programs are initially installed onto the 
hard disk of the authorisation subsystem from CD-ROM. 

35 An input/output interface (I/O) 205 controls interaction of a space 

controller with the authorisation subsystem. A network interface (NET) 210 
processes a set of protocol layers, which allows the authorisation 
subsystem to exchange information with the reservation subsystem (through 
the communication channel) . 

40 



The working memory 155 stores a historical repository 215 for each 
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flight. The historical repository 215 logs information about a series of 
previous instances of the flight; preferably, the last 104 previous 
instances of the flight (one per week in two years) are logged. 
Particularly, the historical repository 215 includes a field INFOj (with 
5 j=l~104) that identifies the previous instance of the flight (for example 
its number and departure time) , and a field SPR, with information about 
each shipment (for example the weight and the volume of the freights that 
have been transported, and a fare actually applied) . The historical 
repository 215 further stores a set of historical profiles, which consist 

10 of a temporal profile of the weight and of the volume that have been 
reserved by each category during an opening period of the previous 
instance of the flight. Each historical weight profile and each historical 
volume profile is defined by a series of daily snapshots (for example 
taken at the end of the day); if the previous instance of the flight was 

15 open 25 days before its departure, the k-th snapshot (with k=-25~0) for 
the i-th category in the j-th historical weight profile consists of a 
-reserved weight Rp ijk and a residual weight Dp ijk , whereas the sante snapshot 
in the j^th historical volume profile consists of a reserved volume Ryij k 
and a residual volume Dv ijk ; the first snapshot gives the reserved capacity 

20 RPij-25/Rv i:J _25 and the residual capacity Dp i:j _25#I>Vij-25 at the opening of the 
previous instance of the flight, whereas a last snapshot gives the 
reserved capacity Rpijo,RVi j0 and the residual capacity Dp ijo ,DVijo at its 
departure . 

25 Th e working memory 155 further stores a current repository 220 for 

the future instance of each flight. Particularly, the current repository 
220 includes a field INFO that describes the future instance of the flight 
(for example the offered capacity P to t,V tot and its scheduled departure 
time) , a field with the corresponding booking information BK, and a field 

30 with a current time t at which the information is taken. The current 

repository^ 220 further stores a set of current profiles, which consist of 
a temporal profile of the weight and of the volume that have been 
currently reserved by each category for the future instance of the flight. 
Each current weight profile and each current volume profile is defined by 

35 a series of daily snapshots. The k-snapshot (with k ranging from -25 to 
the current time t) for the i-th category in the current weight profile 
consists of a reserved weight Rp iJc and a residual weight Dp ik , whereas the 
same snapshot in the current volume profile consists of a reserved volume 
Rv ik and a residual volume Dv ik ; the last snapshot gives the reserved 

40 capacity Rp it ,Rv it and the residual capacity Dp it ,Dv it at the current time t. 
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The booking information BK, the current reserved capacities RPi t ,Rv it 
and the current residual capacities Pp itf Dv it are received from the 
reservation subsystem through the network interface 210. The set of 
profiles stored in the current repository 220 for each flight are copied 
5 to the corresponding set of profiles in the historical repository 215 
after departure of the flight; the shipment field SHP-j is continually 
updated by the input/output interface 205 as soon as the corresponding 

•i 

information becomes available. 

10 The historical repository 215 and the current repository 220 are 

accessed by a weighting module 225, which assigns a weight Wj to each 
previous instance of the flight. The weight Wj depends on a relevance of 
the previous instance of the flight with respect to its future instance. 

15 The weights Wj are input to an estimation module 230, which further 

« 

accesses the historical repository 215 and the current repository 220. The 
estimation module 230 determines a potential (weight and volume) % profile 
from each « corresponding historical profile. The potential profile is 
obtained adjusting the historical profile according to the respective 
20 current profile and an estimated unconstrained (or untruncated) demand of 
the respective capacity variable for the category in the old instance of 
the flight. The unconstrained demand takes into account the periods of 

r 

time in which the category has been closed to sales (because the reserved 
weight and/or the reserved volume for the category have reached their 
25 booking limits defined by the respective nesting authorisations) , and 

estimates a demand not limited by the offered capacity. The estimation of 
the unconstrained demand is carried out independently for the weight and 
for the volume. Each potential profile consists of 25 snapshots; the k-th 
snapshot for the i-th category in the j-th potential weight profile 

30 consists of a potential weight request Rp ijk , whereas the same snapshot in 
the j-th potential volume profile consists of a potential volume request 
Rvyk . The last snapshot of the potential weight profile and of the 
potential volume profile define a final potential weight request Rp [J0 and 

a final potential volume request Rvyo , respectively, for each i-th 
35 category in the j-th previous instance of the flight. 



The estimation module 230 further determines a potential revenue A 
for each i-th category transported by the j-th previous instance of the 



1 



f 



^20010034GB1 9 

flight. The potential revenue is estimated according to the revenue of the 
category in the previous instance of the flight; this value is further 
adjusted to take into account the revenue of the corresponding capacity 
currently reserved for the future instance of the flight. 

5 

The estimation module 230 outputs a series of historical scenarios 
for each flight. The j-th scenario consists of the corresponding final 

potential weight request Rp ij0 , final potential volume request Rvyo , and 

potential revenue A^ for each i-th category. The space controller defines 
10 a user scenario through the input/output interface 205. The user scenario 
consists of a potential weight request Rp f , a potential volume request 

Rv$ , and a potential revenue A, , which are independently forecast by the 

space- controller for each i-th category. A weight w is assigned to the 
user scenario by the space controller through the input/output interface 

is 205. x % 

? » ■ * 

The weight Wj of each historical scenario {Rp iJ0 , Rvyo , A fJ } and the 

weight w of the user scenario {Rp (t Rvt , A f . }, which will be denoted as a 
whole with W a in the following, are input to a normalisation module 235. 
20 The normalising module 235 determines a probability W s of each (historical 
and user) scenario. 

The probabilities W s are input to a nesting order module 240, which 
also receives the potential revenues Ay of each historical scenario from 

25 the estimation module 230 and the potential revenues A, of the user 

scenario from the input/output interface 205. An input weight nesting 
order IN_ORDp and an input volume nesting order IN_ORDv are provided by 
the space controller thorough the input/output interface 205, and then 
transmitted to the nesting order module 240. The nesting order module 240 
30 outputs the nesting orders ORDp,ORDv for the future instance of the 
flight. 



The nesting orders ORDp,ORDv are provided to an optimisation module 
245; the optimisation module 245 further receives the historical scenarios 
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{Rp ij0 , RVfjo , A fJ } from the estimation module 230 and the user scenario 
{Rp f , Rvi , A f } from the input/output interface 205, together with the 

respective probabilities W s from the normalisation module 235. The space 

controller supplies an aggressiveness parameter AG, indicative of an 
5 attitude to the risk, to the optimisation module 245 through the 

input/output interface 205. The optimisation module 245 employs stochastic- 
programming techniques to determine the (optimal) nesting authorisations 

Api,AVi that maximise a potential revenue A of the future instance of the 
flight. 

10 

Stochastic programming techniques are commonly used as a support in 
decision-making processes under uncertainty, such as in financial trading. 
An example of stochastic programming technique is described in *IBM 
Stochastic Solution", Tito A. Ciriani, Stefano Gliozzi, AIROnews Bulletin 
15 of AIRO, VI, n.l - Spring 2000, pp. 6-10. In particular, the module 245 
optimised ' each scenario using a deterministic linear programming 
technique; the results are then combined weighting each scenario according 
to the respective probability. 

20 Similar considerations apply if a whole application program (defined 

by the different modules described above) and the corresponding data are 
structured in a different manner, if other functions are provided, if the 
historical and current repository are replaced by equivalent memory 
structures, and so on. 

25 

With reference now to Figg.3a-3d, the application program when 
running on the authorisation subsystem performs a method 300. The method 
starts at block 302 and then passes to block 304, wherein a menu with a 
series of possible choices is displayed on the monitor of the 
30 authorisation subsystem. 

The method carries out the operations corresponding to the selected 
choice. Particularly, if the space controller has selected a management 
function, the blocks 305-307 are executed, whereas if the space controller 
35 has selected an optimisation function, the blocks 308-374 are executed; 

conversely, if the space controller has selected an exit option the method 
ends at the final block 376. 

Considering now block 305 (management function) , the information 
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about the future instance of the flight is received from the reservation 
subsystem, and it is stored into the current repository. The method passes 
to block 306, wherein the set of profiles stored in the current repository 
are copied to the historical repository after departure of the flight. 
Continuing to block 307, the information about a selected shipment for the 
previous instance of the flight is updated as soon as it is available. The 
method then returns to block 304 waiting for a new command. 

With reference to block 308 (optimisation function) , a flight is 
selected either automatically (in a batch processing of all the flights of 
the airline) or by the space controller (in an interactive processing) . 

The method then assigns the corresponding weight Wj to the previous 
instances of the flight. Particularly, a season coefficient wsj for each j- 
th previous instance of the flight is calculated at block 308. The 
coefficient ws^ depends on a season of the previous instance of the flight, 
that is. on a difference between a period of the previous instancfe of the 
flight afid its future instance; the more the previous instance of the 
flight has happened in a similar period of the year as its future 
instance, the more the season coefficient wsj is bigger. The season 
coefficient wsj is obtained from a matrix having 52 rows and 52 columns , 
which is preferably symmetric; each element [u,v] of the matrix (with 
u,v=l~52) stores the season coefficient wsj of a previous instance of the 
flight that departed in the u-th week of the? year, for a future instance 
of the flight that will depart in the v-th week of the year; the season 
coefficients ws^ have a value ranging from 0 to 1, with all the diagonal 
elements of the matrix (u=v) that are equal to 1. 

Descending into block 312, a time distance coefficient wg., for each 
j-th previous instance of the flight is calculated. The coefficient wg^ 
depends on^a space between the current time and the departure of the 
previous instance of the flight; the more recent the previous instance of 
the flight is, the bigger the time distance coefficient wg.j is. The time 
distance coefficient wgj is defined as a hyperbolic function of a number of 
days g that have passed from the departure of the previous instance of the 
flight. We set wg j= l for g=0 and wgj=a x for g-x* (with 0<a x <l, and 
preferably a^O) ; moreover, we set an exponent a 2 of the exponential law 
(with a 2 <l, and preferably a 2 =0.2), and a scale coefficient a 3 for reducing 
a concavity of the hyperbola (with a 3 >30, and preferably a 3 =30) . The time 
distance coefficient wgj is then given by: 
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This results in a time distance coefficient wgjSO.53 for g=708 (2 
years) , so as to give more relevance to the most recent previous instances 
of the flight, while maintaining some* relevance for the oldest previous 
instances of the flight at the same time. 

The weight w-, for each j-th previous instance of the flight is then 
calculated at block 314, by combining the corresponding season coefficient 
wsj and time distance coefficient wgj . Particularly, the value w^ is 
calculated as a sum of the coefficients wsj and wg j# which are weighted 
according to a parameter h with a value between 0 and 1, and preferably 
equal to 0.65: 

w/= h • wSj +(l-/r)-wg y 

Considering now block 316, an opening coefficient a ijk for each k-th 
snapshot of the i-th category in the j-th previous instance of the flight 
is estimated; the opening coefficient a ijk indicates a percentage of time 
during which the category has been open to sales during a period 
associated with the snapshot, that is from a preceding snapshot (k-1) to a 
current snapshot (k) . 



As shown in Fig. 4a, the category is open only when the residual 
weight and the residual volume are both positive. We denoted with s*(f) a 
function that determines the sign of its argument f, that is: 



li<y/>o 



The starting conditions that are used (in the given order) to 
determine whether the category has been always open or always closed are: 



a ijk = < 



1 '/ s*(t>p yk _ l ) = s + (Dp 0t ) = s+iDv^) = , + (Dv, t ) = 1 

0 '/ * + 0>/W = s\Dp IJk ) = 0 or s+iDv^) = s+(Dv gk ) = 0 
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Otherwise, we calculate the sign of the derivative of the functions 
defined by the residual weight and the residual volume: 



sp' iJk =S^Dp gk )-s + (Dp gl _ l ) 
sv' sk =s+(Dv sk )-s*(Dv ffk _ l ) 



When sp iJk ,SVy k & the intersection of these functions with a time 
axis is given by: 

DPyk-l-DPijk 

^ak-i 



9 " Dv ijk . x -Dv iJk 



We now define the following functions: 



IPijk = 



*Pyi if Wyk =-1 
l-tv (Jk if sv\ jk = +1 



If we consider all the possible cases, we can demonstrate that: 



tPgk if s Pijk ^Oandsvp = 0 
K if ^ijk * 0 and sp] Jk = 0 
minfc^ }// sp] jk -sv' yk > 0 
^ max^),(^ +tv iJk )-l}otherwise 

The method continues to block 318, wherein a weight emphasis value 
Epik and a volume emphasis value Ev ik to be applied to the k-th snapshot for 
each i-th category (irrespective of the old instance of the flight) are 
calculated; the emphasis values Ep ik ,Ev i)c indicate what could have happened 
if the category had been always open to sales in the period associated 
with the k-th snapshot. The emphasis is calculated as a function of 
observations consisting of the demand of the respective category in the 
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corresponding historical profiles where the category was open in the same 
period. 

The method takes into account all the corresponding (weight or 
5 volume) historical profiles for the same category, which were always open 
in the period (a ii)c =l) . Further historical profiles (with a decreasing value 
of the opening coefficient aij k ) are also taken into account, if necessary, 
in order to ensure a sufficient number of observations (for example at 
least equal to 30% of the number of the previous instances of the flight) . 

10 If we denote with Pi k = P gk — P gk ^ and with V- Jk — V gk — V gk _ x a gradient of the 

reserved weight and of the reserved volume, respectively, the emphasis 
values Epik, Ev iIc are defined as a weighted mean of these values: 



15 




/ 

20 

r 

Passing to block 320, the method estimates a potential weight 
gradient Upij* and a potential volume gradient Uvij* for each k-th snapshot 
of the i-th category in the j-th historical weight profile and historical 
volume profile, respectively. The potential gradients Upij*, Uvi-jk are 

25 obtained as a linear interpolation between the value Pgk 9 V ijk for a^sl and 
the value max[f|^,^ A ]max[F^,£v^] for a ijk =0, that is: 

Up gk = m*x[Py k , Ep iJk + a ijk + (P gk - Ep gk )] 
Uv ijk = max[^ , Ev ijk + a yk + (V!., - Ev ijk )] 



The method then builds each potential profile. As shown in Fig: 4b, 
the potential profile is set to the respective current profile up to a 
35 time t' corresponding to the current time t. Considering now block 322, 
the next snapshot for the i-th category in the j-th potential weight 



30 
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profile is calculated integrating the potential weight gradient Up ijJc , 
whereas the same snapshot for the i-th category in the j-th volume profile 
is calculated integrating the potential volume gradient Uv ij)c ; the 

integration is carried out with a simple sum, that is Rp ijk = Rp ijk _ x +Up ijk 

and Rvgt =/?VyjM +Uv ijk , respectively. The potential gradients relating to 
the period including the time t' are reduced to take into account the time 
that has passed from the preceding snapshot; we denote with k an index of 
the snapshot immediately following the time t'; if g t and gj are the 

number of hours from the time t' and from the £-th snapshot to the 
departure of the previous instance of the flight, respectively, the 
potential gradients Up y - 9 Uv g - are reduced according to the following 
formulas : 

up-^up g *- g < 



The method then checks at block 324 whether either the resulting 
potential weight request P gk or the resulting potential volume request V^ k 
for the current snapshot is not strictly positive (>0) . If so, the method 
descends into block 326, wherein the potential weight request P ijk and the 

potential volume request V ijk are both set to zero; the method then 
continues to block 328. Conversely, if both the resulting potential weight 
request P ijk and the resulting potential volume request V ijk are strictly 
positive, the method passes to block 328 directly. In other words, the 
potential weight capacity P ijk and the potential volume capacity V ijk (with 
k=t,&,~0) for the current snapshot are given by: 



Pyk-i + Up IJk if P ijk -\ + Up iJk > 0 and V (J k-\ + Uv gt > 0 
0 otherwise 

V,^ + Uv yk if Pyk-t + Up ijk > 0 and V yk -i + Uv yk > 0 
0 otherwise 
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With reference now to block 328, the result is further corrected by 
reconciling the potential weight request Py k and the potential volume 

request Vg k of the current snapshot, in order to resolve any inconsistency 

that may be introduced between the two capacity variables and that may 
5 bring about incongruent values of a density of the potential request. We 

denote with = the density of the reserved capacity for the i-th 

category in the k-th snapshot of the j-th previous instance of the flight, 

and with Oyk = — — the density of the potential capacity for the same 

snapshot. If the density 8ijk for the current snapshot is comprised between 

10 5ij k and 8gk-\ , not reconciliation is carried out and the method descends 
into block 330 (described in the following) . 

Otherwise, the potential demand is correct to an admissible reference 
15 density R8 ijk at block 332. The reference density R5 ijk is set to the value 

between 5ij k and Syk-i that is closer to 8yk . Particularly, the potential 



weight request Rp^ k and the potential volume request Rp^ k are reconciled 
according to the following formula: 

ht^^(Py*+VskRS iJk ) 

20 - 

Considering now block 330, the method checks whether the last 
snapshot (associated with the departure of the old instance of the flight) 
25 has been reached. If not, the method returns to block 324 for processing a 
next snapshot. 

Conversely, the method descends into block 333, wherein a historical 
revenue for the i-th category transported by the j-th previous instance 
30 of the flight is calculated by summing a value of the respective 

shipments. The historical revenue Aij may be calculated only when the value 
of all the shipments for the category are known; otherwise, if the value 
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of one or more shipments is not available, the historical revenue A Li for 
the whole category is deemed unknown. In a similar manner, a current 
revenue A L for the i-th category in the future instance of the flight is 
calculated by summing a value of the respective requests accepted at the 
current time t. 

In the cargo flights, the freight fares are defined according to a 
so-called taxable weight P ; the taxable weight P* is a virtual weight of 
the freights that is applied to calculate their value for a given unit 
weight fare. We denote with A* a reference density (for example A*=167 
kg/m 3 ) . If a density 5 of the freights is equal to or higher than A* (heavy 
freights) the taxable weight P # of the freights is equal to their real 
weight; otherwise, if the density 8 of the freights is lower than A* (light 
freights), the taxable weight P* of the freights is obtained multiplying 
their 'volume by the reference density A*: 

P-=max(P,VA*) 

A (weight) unit revenue X is then calculated as A = — 

P* 

We define a historical taxable weight, for each i-th category in the 
j-th previous instance of the flight as Rp* 0 = max(Rp &0 , Rv iJ0 • A* ) , and a 
current taxable weight for each i-th category as Rp[ = max(/2/7 ^ , Rv u • A* ) ; it 
is then possible to calculate a historical unit revenue X Lj for each i-th 
category in the j-th previous instance of the flight and a current unit 
revenue A* for each i-th category in the future instance of the flight by 

applying the above mentioned formula, that is A. =-^- and A. =-^- 

respectively. 

Proceeding to block 334, the historical unit revenues X ijt for each 
i-th category in the j-th previous instance of the flight wherein the 
respective historical revenue is unknown, are set to the weighted mean 
of the respective values X i:j available, that is: 



A„. = -A== 



18 



The method then estimates a potential unit revenue Ay for each i-th 
category in the j-th previous instance of the flight, by adjusting the 
corresponding historical unit revenue according to the corresponding 
current unit revenue A*. Particularly, a departure coefficient wd for the 
5 future instance of the flight is calculated at block 336. The coefficient 
wd depends on a space between the current time t and the departure of the 
future instance of the flight. The more the current time is close to the - 
departure, the bigger the coefficient wd is. We assume that the current 
unit revenue X± does not affect the historical unit revenue X ±i when the 
10 number of days to the departure of the future instance of the flight is 
greater than a pre-set value g m (for example g m =10) . If we denote with g d 
the number of days from the current time t to the departure of the future 
instance of the flight, the departure coefficient wd is given by: 




Passing to block 338, an increment coefficient wx^ for each i-th 
category in the j-th previous instance of the flight is further 
20 calculated. The coefficient wxij depends on an increment of the potential 
capacity at the departure for the i-th category in the j-th previous 
instance of the flight with respect to the current reserved capacity. We 
define a potential taxable weight for each i-th category in the j-th 

previous instance of the flight as P^ 0 =max{Rp.j 09 Rp iJ0 - A ); the increment 
25 coefficient wxij is then given by the following formula: 



30 A corrective factor wcij for each i-th category in the j-th previous 

instance of the flight is then calculated at block 340, by combining the 
corresponding departure coefficient wd and increment coefficient wx^. 
Particularly, the corrective factor wcj is calculated as a sum of the 
coefficients wd and wxij, which are weighted according to a parameter he 

35 having a value between 0.5 and 1, and preferably equal to 0.6: 
wc ij =hc wx ij + (1-hc) wd 



max 



/ — r . \ 
P -P 



0 otherwise 
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The potential unit revenue X i} for the i-th category in the j-th old 
instance of the flight is then calculated at block 342 as a sum of the 
corresponding historical unit revenue and current unit revenue which 
are weighted according to the corrective factor wc^, that is: 



The potential revenue A^. for each i-th category in the j-th old 
instance of the flight is simply obtained as = - Rp* g0 . 

Proceeding to block 344, the historical scenarios { Rp ij0 , Rvijo , XT } 
are determined for each old instance of the flight. The method then checks 
at block 348 whether the space controller wishes to input a user scenario. 
If not, 'the method descends into block 350 (described in the following) . 

Conversely, the space controller inputs the user scenario {Rp t ,Rv ( , } at 

block 352 and inputs its weight w at block 354; the method then passes to 
block 350. 



Considering now block 350, the probability w s of each (historical 

and user) scenario is calculated normalising the corresponding weight w s , 
that is: 

W. = 1 — 



* \04jver 

I 



The method then proceeds to block 355, wherein the input nesting 
orders IN_0RDp, IN.ORDv are provided by the space controller. Descending 
into- block 356, the method checks whether for each i-th category some of 
the j-th historical scenarios have the final potential weight request 

RPgo ' the final potential volume request Rvyo and the potential revenue 

all strictly positive. If not, the method passes to bock 358 
(described in the following). Conversely, the method enters block 360, 



iil*i0010034GBl 



20 



15 



wherein a potential weight revenue = and a potential volume 

revenue AV< J =-=^^ are calculated for each i-th category that meets the 
above mentioned conditions (taking into consideration all the historical 

A 

scenarios and the user scenario) . A corresponding mean weight revenue Ap f 

A 

and a corresponding mean volume revenue Xv { are then calculated for these 
categories as a weighted mean of the potential weight revenues *kp^ and the 
potential volume revenues Xv b , respectively: 



10 ^=-h== 

j 

s 



The method then passes to block 358. 



Considering now block 358, the position of each category in the 
nesting orders ORDp,ORDv is then determined. For each capacity variable 
20 (weight and volume) , the categories that do not have a respective mean 
revenue maintain their position in the input nesting orders 
IN_ORDp, IN_ORDv. The other categories having a respective mean revenue 

A A 

atife ranked in a corresponding descending order within the 
available positions. For example, let us suppose that the input weight 
25 nesting order for the five categories is IN_ORDp={C 2 ,C4,C3,C 5 ,Ci}. The 
categories C 2 \ C 3 and C 5 are ranked in the order {C 3 ,C 5 ,C 2 } according to 

AAA 

their respective mean weight revenues Ap 3 , Ap $ , Ap x ; conversely, the mean 

weight revenue Xp A (for the category C 4 ) and the mean weight revenue 7<p x 
(for the category Ci) are not available. In this case, the resulting weight 
30 nesting order is ORDp={C 3 , C 4 ,C 5 ,C 2 ,C 1 }. 
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The method then proceeds to block 362, wherein the aggressiveness 
parameter AG is input by the space controller. The method builds a 
stochastic model that maximises a target function defined by the potential 

revenue A of the future instance of the flight; the potential revenue A 
is calculated using both a nesting policy (irrespective of an order of the 
requests) and an out of nesting policy; the nesting policy (best case) and 
the out of nesting policy (worst case) are combined according to the 
aggressiveness parameter. 

Particularly, a portion of the target function for the nesting policy 
is built at block 364. We define a potential nesting constrained weight 

TNpy and a potential nesting constrained volume TNVy that would have been 

transported for each i-th category by the j-th previous instance of the 
flight complying with the corresponding nesting authorisations Ap ± and Av i# 
respectively: 

^mr 9 <Av ( 



The variables TNp~ , TNv.j are subjected to the constraint 

► 77VVT > ► 

TNp g = ^=r- and to the simple bounds TNp g < P off , TNv g < V off . 

S ijO 



A potential nesting revenue TNAj of the j-th previous instance of 
the flight that would have been produced by the potential nesting 

constrain^ weight TNp (j and the potential nesting constrained volume TNv^ 
is given by: 

5 , JL^ 

mA J = y £mp ij .* Pi 

The method then builds a further portion of the target function for 
the out of nesting policy, using two independent components for the weight 
and for the volume. Particularly, a set of out of nesting weight 
authorisations AOpi and a set of out of nesting volume authorisations AOVi 
are defined at block 366. If we denote with Pe the offered weight that is 
not used, and with Ve the offered volume that is not used, the out of 
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nesting authorisations AOpi,AOVi are defined by the following relations: 

5 

P = 1 £AOp i +Pe 

r=l 

V^AChi+Ve 

1=1 



The out of nesting authorisations AOpi and AOVi are further 
constrained by the nesting authorisations Ap A and Av i# respectively: 

5 

J^Apt^AOpt+Pe 

i=\ & 

j^Av^J^AOVt+Ve 

i=l 4 * 

We now define a potential out of nesting constrained weight TOpg 

that would have been transported for each i-th category by the j-th 
previous instance of the flight in the worst case complying with the out 
of nesting authorisations AOpi,AOVi only for the weight, that is: 



TOp^AOp, 
fd^<AOv r X^ 

Even in this case, the potential out of nesting constrained weight 
TOpg is subjected to the simple bounds TOpg £ P off . A potential out of 

nesting weight revenue TOAp. of the j-th previous instance of the- flight 
that would have been produced by the potential out of nesting constrained 

weight Topy is given by: 

5 , A 

TOApj^TOpyXp, 

Proceeding to block 368, we now define a potential out of nesting 

constrained volume TOVg that would have been transported for each i-th 
category by the j-th previous instance of the flight in the worst case 



0010034GB1 



23 



complying with the out of nesting authorisations AOpi.AOvj. only for the 
volume, that is: 



TOv g < AOv. 

As above, the potential out of nesting constrained volume TOVg is 
subjected to the simple bounds TOv i} < V M . A potential out of nesting 

volume revenue TOAVj of the j-th previous instance of the flight that 
would have been produced by the potential out of nesting constrained 
volume Tov {J is given by: 

The method then passes to block 370, wherein the target function to 
be optimised is defined combining the nesting portion and the weight and 
volume component of the out of nesting portion with the respective 
probabilities, according to the aggressiveness parameter AG: 

3 S 2 * 

The nesting authorisations Ap^Avj. that maximise the target function 
defined above (with the variables meeting the respecting constrains and 
single bounds) are then obtained at block 372. 

Several by-products of this process may also be obtained. For 
example, the stochastic model calculates the minimum selling price MXp,MXv; 
particularly, the value MAp is equal to a marginal weight revenue and the 
value MXv is equal to a marginal volume, revenue in the neighbourhoods of 
the optimal solution, which are obtained as the partial derivative of the 
target function with respect to the offered weight and offered volume, 
respectively. Moreover, the stochastic model calculates a mean potential 

constrained weight CRp g and a mean potential constrained volume CRv g as: 
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AG^(TNp s w s ) + (l-AG)^(X0p s w s ) 

S S 

AG^(7^sV t )+(l-AG)^(fd7,^,) 

s s 

The method verifies at block 374 whether a last flight to be 
optimised has been processed. If not, the method returns to block 308 for 
applying the operations described above to a next flight. Conversely, the 
method returns to block 304 waiting for a new command. 

Similar considerations apply if the program application performs an 
equivalent method, for example with error routines, escape functions for 
stopping execution of the program, several concurrent processes that . 
execute the above described operations in parallel, and the like. 
Alternatively, a resolution of the freights or a unit load device used for 
transporting the freights are considered (for example by reducing the 
residual ^capacity of a pre-set coefficient), no-show freights and/or 
denied boarding (off-loaded) freights are estimated, an overbooking is 
taken into account (for example increasing the offered capacity by a pre- 
set coefficient), or no by-product is calculated. 

More generally, the present invention provides a yield management 
method, which is used to optimise a yield parameter resulting from 
assigning a capacity offered by a future instance of a service to each one 
of a plurality of different categories of requests competing for the 
capacity; the capacity is defined by one or more capacity variables. A set 
of historical profiles is stored for each one of a plurality of previous 
instances of the service; the set includes a historical profile of an 
historical value of each capacity variable reserved by each category. A 
probability is then assigned to each previous instance of the service. The 
method estimates a potential profile of a potential value of the capacity 
variable from each historical profile; the estimation is carried out 
according to a corresponding current value of the capacity variable 
reserved by the category for the future instance of the service and 
according to a corresponding unconstrained demand of the capacity variable 
for the category in the previous instance of the service. A historical 
scenario is defined for each previous instance of the service; the 
historical scenario includes a final potential capacity variable from each 
corresponding potential profile. The method then determines an 
authorisation to allocate the offered capacity for each capacity variable 



CRPy = 
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of each category in the future instance of the service by applying a 
stochastic model to the historical scenarios according to the 
corresponding probabilities. 

5 The computational complexity of the method according to the invention 

increases linearly with the number of the historical scenarios and it is 
independent of the number of categories; this is particularly advantageous 
in a commercial application wherein the market dynamic is very fast, so 
that a relatively low number of historical scenarios is often enough to 
10 take planning decisions. As a consequence, the problem defined by the 

yield management method can be solved in a reasonable amount of time. This 
allows the nesting authorisations to be recalculated several times during 
the opening period, and particularly with a higher frequency close to the 
departure of the future instance of the flight. 

15 

The solution of the invention does not require any simplification, 
shortcut or assumption to be made. Particularly, all the capacitiy 
variable^ are taken into account simultaneously, the scenarios are not 
artificially reduced in any way, and no assumption is made on the (simple 
20 and joint) statistical distribution of the variables involved in the 
problem. , 

The envisaged solution derives the probabilistic nature of the demand 
directly from the reality. As a consequence, any kind of demand (even if 
25 particularly scattered or inconsistent) yields to consistent and 

conservative suggestions. Therefore, the method of the invention is 
accurate, robust and secure. 

The preferred embodiment of the invention described above offers 
30 further advantages. For example, the probability assigned to each 

historical scenario is calculated by normalising the corresponding weight, 
which is a combination of the season coefficient and the time distance 
coefficient. In this way, the previous instances of the flight that have 
happened in a similar period as their future instance and/or that are 
35 close to the^ current time are given more importance. 

Preferably, this value is calculated as a weighted sum of the season 
coefficient and the time distance coefficient. This approach gives more 
relevance to the dominant coefficient, and yields to a good result even if 
40 the number of previous instances of the flight is not very high. 
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Similar considerations apply if a different number of periods is 
taken into considerations (such as the days of the week) , if the season 
coefficient and/or the time distance coefficient are calculated in a 
different manner, if different values are used in the formula defining the 
5 time distance coefficient and/or the weight, if the time distance 

coefficient depends on a release of the previous instance of the flight 
(or more generally on an occurrence thereof) , and so on. Alternatively, 
each weight is assigned to a set consisting of more than, one previous 
instance of the flight (such as including all the previous instances of 
10 the flight in a corresponding month) , the weight further depends on mobile 
events (such as Easter), the weight is calculated as the product of the 
season coefficient and the time distance coefficient (giving a high weight 
only when both the coefficients are high) . 

15 The potential profiles are estimated by emphasising the demand when 

the respective category has been closed to sales; the emphasis values are 
defined in order to take into account both the weight and the volume in a 
symmetric. manner. The suggested procedure emphasises the demand according 
to what has happened in the other previous instances of the flight in 

20 which the category was (almost) always open. 

Moreover, the result of the integration is corrected to ensure that 
no negative forecasts for the potential requests at the departure of the 
future instance of the flight are obtained when the demand is low. 

25 

Similar considerations apply if each profile consists of a different 
number of snapshots, if the snapshots are taken with a different 
frequency, if the potential profiles are built with a different procedure, 
if the opening coefficients and/or the emphasis values are' calculated in 
30 another manner, if a different interpolation technique for obtaining the 

potential gradients is employed, if the potential gradients are integrated 
in a different manner, if a dummy snapshot preceding the opening of the 
flight is introduced for processing flights that have not been opened yet, 
and the like. . 

35 

Alternatively, the emphasis values are calculated taking into account 
only the previous instances of the flight in which the category was 
always open (assuming a sufficient number of corresponding historical 
profiles) , using only the historical profiles with a growing demand, or 
40 performing a weighted mean of the more optimistic increments of the demand 
in the historical profiles. However, the method according to the present 
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invention leads itself to be carried out without any correction to the 
result of the integration, or even detecting and storing any denied 
request, so as to know the potential profiles directly (even if some 
adjustment is always suggested because it is very difficult to known all 
the denied requests for sure) . 

The potential profiles are further reconciled to an admissible 
reference density; in this way, it is assured that no inconsistent 
forecasts are produced. 

The reference density is set to the closest one between the density 
of the previous snapshot in the potential profile and the density of the 
current snapshot in the historical profile.. This value makes it possible 
to take into account what has happened in the historical profile without 
excessively affecting the potential profile. 

Similar considerations apply if another logic relation between the 
capacity^ variables is considered for the reconciliation of the potential 
profiles, if the reconciliation is carried out in a different manner, if a 
different reference density is employed, and so on. For example, in an 
alternative embodiment of the present invention the reference density is 
set to a mean of the above-mentioned limit values, to the density of the 
preceding snapshot in the historical profile, or to a pre-set value. 
However, the method of the invention leads itself to be carried out even 
without performing any reconciliation of the potential profiles. 

Each historical scenario further includes the potential revenue, 
which is estimated from the corresponding historical revenue; when the 
value is unknown, the historical revenue is set to a weighted mean of the 
values known in the other previous instances of the flight. This is 
particularly advantageous when it is difficult, if not impossible, to 
detect trie real revenue of each shipment in the previous instance of the 
flight at its departure; moreover, the employed procedure is independent 
of the hypothesis that the freights of each shipment maintain the same 
category from the historical profile to the potential profile. 

Preferably, the potential unit revenue is calculated as a mean of the 
historical unit revenue and the current unit revenue, weighted according 
to the corrective factor. In this way, the potential unit revenue takes 
into account what is currently happening for the future instance of the 
flight; however, the current unit revenue affects the potential unit 
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revenue only when the corrective factor is high (this is preferable when 
the current revenue is not particularly accurate) . 

The corrective factor is a combination of the departure coefficient 
5 and the increment coefficient. In this way, the potential revenue takes 
into account what is actually happening only when this can affect the 
forecast . 

The corrective factor is calculated as a weighted sum of the 
10 departure coefficient and the increment coefficient. This approach gives 
more relevance to the dominant coefficient, and yields to a good result 
even if the current revenue is not very accurate. 

Similar considerations apply if the potential revenue is estimated in 
15 a different manner, if the value of each accepted request is set to a mean 
market value (when its value cannot be determined in an accurate manner by 
the reservation subsystem), if another formula is employed for C9mbining 
the potential unit revenue and the current unit revenue, if the departure 
coefficient and/or the increment coefficient are calculated in a different 
20 manner, if different values are used in the formula defining the departure 
coefficient and/or the corrective factor, and so on. Alternatively, the 
corrective factor depends on other coefficients, the corrective factor is 
calculated as the product of the departure coefficient and the increment 
coefficient (giving a high value only when both the coefficients are 
25 high) ; however, the solution of the invention leads itself to be 

implemented without taking the current revenue into consideration, or even 
without estimating any potential revenue. 

Two nesting orders (for the weight and the volume) aire calculated 
30 according to the mean potential revenue of each category. 

Moreover, corresponding input nesting orders are also provided; the 
nesting orders are obtained updating the input nesting orders according to 
the mean potential revenue for the categories in which this value is 
35 available. 

Similar considerations apply if another technique is employed for 
calculating the nesting orders, or if the input nesting order are updated 
in a different manner. Alternatively, no input nesting orders are 
40 provided, or the nesting orders are not calculated by the system (but they 
are always set by the space controller) . 
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The nesting portion and the out of nesting portion of the target 
function, which are combined using the aggressiveness parameter, allow the 
result of the optimisation to be tuned according to different attitudes to 
the risk; particularly, a more aggressive behaviour gives a higher revenue 
in the worst case (wherein the requests for the most valuable freights are 
the last ones) , and a lower revenue in the best case (wherein the requests 
for the most valuable freights are the first ones) . 

Advantageously, the out of nesting portion of the target function is 
calculated with two independent components for the weight and for the 
volume. In this way, the constraints for the worst case are slightly 
relaxed with respect to a real out of nesting policy, wherein the out of 
nesting authorisations must be met by both the weight and the volume. 
Conversely, in the proposed formula the two capacity variables are 
considered one at the time, so that each request may be accepted even if 
it meets the out of nesting authorisations only for the weight and not for 
the volume (or vice-versa) : % 

Similar considerations apply if the stochastic model is defined in a 
different manner, if another target function is envisaged, if the 
different .components of the out of nesting portion and/or the two portions 
of the target function are combined in another way, and the like. However, 

r 

the solution of the invention leads itself to be carried out even with the 
out of nesting portion that is calculated simultaneously for the weight 
and the volume, or with the target function that is defined with a single 
policy (either in nesting or out of nesting) . 

Preferably, the space controller inputs a user scenario with the 
respective probability; in this way, any personal feeling of the future 
may be introduced into the stochastic model in a very simple manner. 
However, the solution according to the present invention leads itself to 
be implemented even without the possibility of imputing any further 
scenario. 

The proposed solution is particularly advantageous in a cargo flight, 
wherein the capacity is defined by two independent variables (weight and 
volume) , and the revenue of the old instances of the flight is not 
immediately available. Alternatively, the solution of the invention is 
used to optimise another yield parameter resulting from assigning an 
offered capacity defined by one or more different variables; moreover, 
other applications of the yield management method are contemplated, for 
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example in a passenger flight, in a different transportation system (such 
as a train), in a system for controlling booking of rooms in a hotel, and 
the like. 

Advantageously, the solution according to the present invention is 
implemented with a program application, which consists of a computer 
program (software) provided on CD-ROM. 

Alternatively, the program application is provided on floppy-disk, is 
pre-loaded onto the hard-disk, or is stored on any other computer readable 
medium, is sent to the authorisation subsystem through the INTERNET, is 
broadcast, or more generally is provided in any other form directly 
loadable into the working memory of the authorisation subsystem. However, 
the method according to the present invention leads itself to be carried 
out even with a hardware structure, for example integrated in a chip of 
semiconductor material. 

Naturally, in order to satisfy local and specific requirements, a 
person skilled in the art may apply to the solution described above many 
modifications and alterations all of which, however, are included within 
the scope of protection of the invention as defined by the following 
claims . 
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CLAIMS 

1. A yield management method (300) for optimising a yield parameter 
resulting from assigning a capacity offered by a future instance of a 
service to each one of a plurality of different categories of requests 
competing for the capacity, the capacity being defined by at least one 
capacity variable, the method including the steps of: 

storing (306) a set of historical profiles for each one of a 
plurality of previous instances of the service, the set including a 
historical profile of an historical value of each capacity variable 
reserved by each category, 

assigning (310-314,350) a probability to each previous instance of 
the service, 

estimating (316-332) a potential profile of a potential value of the 
capacity variable from each historical profile according to a 
corresponding current value of the capacity variable reserved by the 
category for the future instance of the service and according tc* a 
corresponding unconstrained demand of the capacity variable for the 
category in the previous instance of the service, 

defining (344) a historical scenario for each previous instance of 
the service, the historical scenario including a final potential capacity 
variable from each corresponding potential profile, and 

determining (362-372) an authorisation to allocate the offered 
capacity for each capacity variable of each category in the future 
instance of the service by applying a stochastic model to the historical 
scenarios according to the corresponding probabilities, 

2. The yield management method (300) according to claim 1, wherein the 
step of assigning the probability (310-314,350) includes: 

determining (310) a first coefficient depending on a difference 
between a temporal period associated with the previous instance of the 
service and a temporal period associated with the future instance of the 
service, 

determining (312) a second coefficient depending on a space between a 
current time* and an occurrence time of the previous instance of the 
service, 

calculating (314) a weight by combining the first coefficient and the 
second coefficient, 

calculating (350) the probability by normalising the weight. 

3. The yield management method (300) according to claim 2, wherein the 
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step of calculating the weight consists of calculating (314) a weighted 
sum of the first coefficient and the second coefficient. 

4. The yield management method (300) according to any claim from 1 to 3, 
5 wherein each historical profile and each potential profile include a 

plurality of corresponding snapshots of the reserved capacity variable and 
of the potential capacity variable, respectively, the step of estimating 
the potential prof ile- (316-332) including: 

estimating (316) an opening coefficient for each period comprised 
10 between two consecutive snapshots, the opening coefficient being 

indicative of a percentage of time during which the category was open in 
the period, 

calculating (318) an emphasis value for each period as a weighted 
mean of a gradient in the period of the reserved capacity variable for the 
15 category in a subset of the corresponding historical profiles, 

estimating (320) a potential gradient for each period as a linear 
interpolation between the gradient, for a first value of the opening 
coefficient indicative of a complete opening of the category, and the 
highest between the gradient and the emphasis value, for a second value of 
20 the opening coefficient indicative of a complete closure of the category, 
and 

constructing (322) the potential profile from a time corresponding . to 
the current time by. integrating the potential gradients starting from the 
corresponding current capacity variable. 

25 

5. The yield management method (300) according to claim 4, wherein the 
step of estimating the potential profile (316-332) further includes: 

verifying (324) whether at least one result of the integration in 
each snapshot for each category is not strictly positive, and 
30 setting (326) each potential capacity variable of the category in the 

snapshot to zero if the verification is affirmative. 

6. The yield management method (300) according to claim 4 or 5, wherein 
the at least one capacity variable consists of a plurality of capacity 

35 variables, the step of estimating the potential profile (316-332) further 
including reconciling (328,332) the potential capacity variables for each 
category in the snapshot to a reference value of a logic relation 
therebetween . 



40 



7. The yield management method (300) according to claim 6, wherein the 
step of reconciling (328,332) includes: 
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verifying (328) whether the logic relation between the potential 
capacity variables for the category in the snapshot is included between a 
first limit defined by the logic relation between the reserved capacity 
variables for the category in the snapshot and a second limit defined by 
the logic relation between the potential capacity variables for the 
category in a preceding snapshot, and 

updating (332) the potential capacity variables for the category in 
the snapshot for correcting the corresponding logic relation to the 
closest one of the first limit and the second limit. 

8. The yield management method (300) according to any claim from 1 to 
claim 7, further including the steps of: 

determining (333) a historical unit value of the yield parameter for 
each category in each previous instance of the service if available, 

estimating (334) the historical unit yield parameter for each 
category in the other previous instances of the service as a weighted mean 
of the corresponding historical unit yield parameters available* 

estimating (305,336-342) a potential unit value of the yield 
parameter for each category in each previous instance of the flight from 
the corresponding historical unit yield parameter, and 

calculating (344) a potential value of the yield parameter for each 
category in each previous instance of the service multiplying the 
corresponding potential unit yield parameter by the corresponding 
potential capacity variable, the potential yield parameter being included 
in the corresponding historical scenario. 

9. The yield management method (300) according to claim 8, wherein the 
step of estimating (305,336-342) each potential unit yield parameter 
includes: 

determining (305) a current unit value of the yield parameter for the 
corresponding category in the future instance of the flight, and 

calculating (336-342) each potential unit yield parameter as a sum of 
the corresponding historical unit yield parameter and current unit yield 
parameter weighted according to a corrective factor. 

10. The yield management method (300) according to claim 9, wherein the 
step of calculating (336-342) the potential unit yield parameter further 
includes: 

determining (336) a further first coefficient depending on a 
difference between the current time and a planned occurrence time of the 
future instance of the service, 
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determining (338) a further second coefficient depending on an 
increment of the at least one potential capacity variable with respect to 
the at least one current capacity variable for the category, 

calculating (340) the corrective factor by combining the further 
first coefficient and the further second coefficient, and 

calculating (342) the potential yield parameter as a sum of the 
historical yield parameter and the current yield parameter weighted 
according to the corrective factor. 



10 11. The yield management method (300) according to any claim from 8 to 
10, further including the steps of: 

calculating (356,360) a weighted mean value of the potential yield 
parameter for each capacity variable of each category, 

determining (355,358) a nesting order of the categories for each 
15 capacity variable according to the corresponding weighted mean potential 
yield parameters. 

12. The. yield management method (300) according to claim 11, wherein the 
step of determining each nesting order (355,358) includes: 

20 providing (355) an input nesting order of the categories for each 

capacity variable, and 

updating (358) each input nesting order by ranking the categories . 
having at lest one scenario with each component thereof that is strictly 
positive. 

25 

13. The yield management method (300) according to any claim from 1 to 
12, wherein the step of determining the authorisations (362-372) includes: 

providing (362) an aggressiveness parameter indicative of an attitude 
to the risk, 

30 defining (364) a first portion of a target function, the first 

portion calculating the yield parameter by assigning the offered capacity 
with a nesting policy, 

defining (366,368) a second portion of the target function, the 
second portion calculating the yield parameter by assigning the offered 
35 capacity with an out of nesting policy, 

defining (370) the target function as a sum of the first portion and 
the second portion weighted according to the aggressiveness parameter, and 

calculating (372) the authorisations by optimising the target 
function. 

40 

14. The yield management method (300) according to claim 13, wherein the 
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step of defining the second portion (366,368) includes defining an 
independent component of the second portion for each capacity variable. 

15. The yield management method (300) according to any claim from 1 to 

14, further including the step (348-354) of providing a user-defined 
scenario with a corresponding probability, the stochastic model being 
further applied to the user-defined scenario according to the 
corresponding probability. . 

16. The yield management method (300) according to any claim from 1 to 

15, wherein the service consists of a cargo flight, the at least one 
capacity variable consists of a weight and a volume, and the yield 
parameter consists of a revenue. 

17.. A computer program application (225-245) directly loadable into a 
working memory (155) of a computer (110) for performing the method of any 
claims from 1 to 16 when the application program is run on the computer. 

18. A program product (170) comprising a computer readable medium on 
which the program application (2225-245) of claim 17 is stored. 

19. A yield management system (100) for optimising a yield parameter 
resulting from assigning a capacity offered by a future instance of a 
service to each one of a plurality of different categories of requests 
competing for the capacity, the capacity being defined by at least one 
capacity variable, the system including means (215) for storing a set of 
historical profiles for each one of a plurality of previous instances of 
the service, the set including a historical profile of an historical value 
of each capacity variable reserved by each category, means (225,235) for 
assigning a probability to each previous instance of the service, means 
(220,230) for estimating a potential profile of a potential value of the 
capacity 'variable from each historical profile according to a 
corresponding current value of the capacity variable reserved by the 
category for the future instance of the service and according to a 
corresponding uncons trained demand of the capacity variable for the 
category in the previous instance of the service, means (230) for defining 
a historical scenario for each previous instance of the service, the 
historical scenario including a final potential capacity variable from 
each corresponding potential profile, and means (240,245) for determining 
an authorisation to allocate the offered capacity for each capacity 
variable of each category in the future instance of the flight by applying 
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a stochastic model to the historical scenarios according to the 
corresponding probabilities. 
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ABSTRACT 

A YIELD MANAGEMENT METHOD AND SYSTEM 

5 A yield management method (300) and system, particularly for 

maximising a revenue that can be obtained from a given capacity that is 
offered by a cargo flight; the capacity is defined by two different 
variables consisting of the weight and the volume of the freights that may 
be transported by the flight. A set of historical profiles of the capacity 

10 (weight and volume) reserved by each category is stored (306) for 

different previous instances of the flight. Corresponding potential 
profiles are estimated (316-326,332) independently for the weight and the 
volume. The estimation is carried out taking into consideration the 
corresponding capacity currently reserved by the category for a future 

15 instance of the flight; moreover, this result is emphasised according to a 
corresponding unconstrained demand for the category (not limited by the 
offered capacity) . A series of historical scenarios of the demand at the 
departure of each previous instance of the flight are then obtained from 
the potential profiles (344). A stochastic model is applied (362-372) to 

20 the historical scenarios according to the corresponding probabilities for 
calculating the optimised weight and volume authorisations for each 
category. 
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